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ABSTRACT 



We present the results of an intense photometric monitoring in the near-infrared (~0.9 yum) with the TRAPPIST robotic telescope 
of the newly discovered binary brown dwarf WISE J104915. 57-531906.1, the third closest system to the Sun at a distance of only 2 
pc. Our twelve nights of photometric time-series reveal a quasi-periodic (P = 4.87 ± O.Olh) variability with a maximal peak-peak 
amplitude of ~ 1 1 % and strong night-to-night evolution. We attribute this variability to the rotational modulation of fast-evolving 
weather patterns in the atmosphere of the coolest component (~Tl-type) of the binary, in agreement with the cloud fragmentation 
mechanism proposed to drive the spectroscopic morphologies of brown dwarfs at the L/T transition. No periodic signal is detected 
for the hottest component (~L8-type). For both brown dwarfs, our data allow us to firmly discard any unique transit during our 
observations for planets > 2R m . For orbital periods smaller than ~9.5h, transiting planets are excluded down to an Earth-size. 
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1. Introduction 

Both inaugurated in 1995 (Rebolo et al. 1995; Mayor & Queloz 
1995), the observational studies of brown dwarfs (BD) and ex- 
oplanets are two of the most active fields of modern astron- 
omy. The atmospheric study of these substellar objects has seen 
tremendous advances , thanks to the sophistication of models 
and the constant improvement of instruments (see, e.g., Seager 
& Deming 2010 and Showman & Kaspi 2012 for reviews). With 
effective temperatures ranging from ~300-2000 K, L and T-types 
field BDs amenable for detailed direct studies represent critical 
precursors to the atmospheric characterization of giant exoplan- 
ets. The data gathered so far outline the important role of at- 
mospheric condensates on the spectral morphologies of these 
objects (Kirckpatrick 2005). This is especially true at the L- 
T transition (~L7-T4 spectral types) which is characterized by 
an increase of the /-band luminosity with decreasing tempera- 
ture (Vrba et al. 2004). This has been explained by the gradual 
depletion of silicates in the cooler atmospheres resulting in in- 
creasingly patchy cloud covers and thus increasingly small con- 
densate opacity. Still, the details of this transition remain poorly 
understood (e.g. Saumon & Marley 2008). Ackerman & Marley 
(2001) have proposed to explain the condensate depletion by the 
fragmentation of the clouds driven by convection in the tropo- 
sphere. This scenario predicts relatively large (~l-20%) photo- 
metric variability around 1/im on rotational timescales, driven by 
the formation, evolution and complex dynamics of cloud holes 
in the upper atmosphere. 

Brown dwarfs being very rapid rotators (typically a few 
hours, Herbst et al. 2007), this variability is observable within a 
few nights of photometric monitoring. Because of the extreme 
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faintness in the optical at the L/T transition, only a handful 
of near-infrared (NIR) monitoring using medium-sized ground- 
based telescopes or space-based facilities could reach the pho- 
tometric precisions required to detect these predicted semi- 
periodic variabilities (Clarke et al. 2008; Artigau et al. 2009; 
Radigan et al. 2012; Khandrika et al. 2013; Apai et al. 2013). 
Some of these results agree nicely with the cloud fragmentation 
hypothesis, but there is still no clear evidence that BDs at the L/T 
transition are more variable than the other BDs (Khandrika et al. 
2013). More high-precision photometric monitoring of L/T tran- 
sition BDs are thus highly desirable. High-precision time-series 
can also inform us on the spatial and temporal distribution of 
cloud structures, vertical thermal profiles, degrees of differen- 
tial rotation, and on the age of brown dwarfs since their rotation 
period decreases monotonically with time. 



A unique opportunity for this domain came up recently with 
the detection by Luhman (2013, hereafter L13) of a nearby bi- 
nary BD at only 2.02 ± 0.15 pc. This system, WISE J104915.57- 
531906.1 (hereafter Luhman 16, following Mamajek 2013), is 
the third closest system to Earth, making it an exquisite target for 
high signal-to-noise ratio follow-up (J = 10.7, K - 8.8). With 
spectral types L8±l for Luhman 16A and Tl±2 for Luhman 16B 
(Kniazev et al. 2013, hereafter K13; see also Burgasser et al. 
2013), both of its components span the L/T transition. This sys- 
tem is an invaluable target for high-precision time-series pho- 
tometry, an interest further reinforced by the classification of 
the pair as a possible variable in WISE All-Sky Source Catalog 
(L13). This motivated us to perform an intensive photometric 
monitoring of the system in the NIR using the 60cm robotic tele- 
scope TRAPPIST. We observed Luhman 16 for twelve nights, 
reaching a precision that allowed us to detect a clear quasi- 
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periodic variability, attributed to the T-dwarf Luhman 16B, and 
combined with a fast evolution of the observed patterns. 

The next section presents our TRAPPIST data and their re- 
duction. Our analysis of the resulting photometric time-series 
is described in Sec. 3. Finally, we discuss briefly our results in 
Sec. 4. 



2. Data description 

We monitored WISE 1049-5319 for twelve nights between 2013 
March 14 and 26 with the robotic 60cm telescope TRAPPIST 
(TKAnsiting Planets and Planetes/mals Small Telescope; Gillon 
et al. 2011, Jehin et al. 2011) located at ESO La Silla 
Observatory in the Atacama Desert, Chile. TRAPPIST is 
equipped with a thermoelectrically-cooled 2K x 2K CCD having 
a pixel scale of 0.65" that translates into a 22' x22' field of view. 
The observations were obtained with an exposure time of 1 15s, 
with the telescope focused and through a special 7 + z' filter 
having a transmittance >90% from 750 nm to beyond 1 100 nrrr 
Considering the transmission curve of this I+z filter, the spectra 
response curve of the CCD detector, and the spectral type of the 
target, we derive an effective wavelength ~910nm for the obser- 
vations. The positions of the stars on the chip were maintained 
to within a few pixels over the course of each run, thanks to a 
'software guiding' system deriving regularly an astrometric so- 
lution for the most recently acquired image and sending pointing 
corrections to the mount if needed. 

After a standard pre-reduction (bias, dark, flatfield correc- 
tion), the stellar fluxes for each run were extracted from the 
images using the IRAF/DAOPHOlj^] aperture photometry soft- 
ware (Stetson, 1987). The same photometric aperture of 8 pix- 
els (~5.1") was used for all nights. After a careful selection of 
stable reference stars of similar brightness, differential photome- 
try was then obtained. Finally, the light curves were normalized. 
The twelve resulting light curves are shown in Fig. 1. To assess 
the night-to-night variability of the target, we also extracted a 
global differential light curve that is shown in Fig. 2. 

The typical full-width at half maximum of TRAPPIST point- 
spread function (PSF) is ~3 pixels = ~2", At this resolution, the 
two components of the 1 .5" binary are only partially resolved, 
and our photometry extracted with an aperture of ~5.1" radius 
shows thus the evolution of the sum of the fluxes of both com- 
ponents. 



^f 






y^f 





366.6 


366.S 


367.6 


367 .B 


1.05 




- 


- 


- 














- r\ 


y*v* 


_ 


_ 


1XJ0 


: 1/ 


V 






0.95 


- 


- 


- 


- 




- i 


... 1 ...... - 


- i , 


i ,,..,, - 






-"\- 



V"' v 





370.6 370.S 


371.6 371. S 




1.05 






- 


1.00 


' / w\ A 

:/ ^ : 




/■ 


0.95 


" 1 1 " 


: 


: 



\J\ r\J 



\ v 



A 



A 



.**v 



V 



/ 



3. Data analysis 

While being relatively stable on longer timescales (Fig. 2), 
Luhman 16 shows a clear variability on a nightly timescale 
(Fig. 1). Furthermore, the observed patterns evolve strongly 
from one night to the other. A Lomb-Scargle periodogram 
(Lomb 1976; Scargle 1982) applied to our photometry shows 
a strong power excess around ~0.2d, matching well the typical 
separation between similar features in the light curves. 

We performed a global analysis of our twelve light curves, 
adapting for that purpose the Markov Chain Monte Carlo 
(MCMC) code described in Gillon et al. (2012, hereafter G12). 
In a first step, we assumed that the observed variability was origi- 
nating from only one of the two BDs. Our basic model for the ro- 



1 http://www.astrodon.com/products/filters/near-infrared/ 

2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 



•••-' 



••. 



>■ 



. 



\ J * ff 



\ / \ 

V 



J#r 



V 



HJD - 2 456 000 



Fig. 1. Normalized TRAPPIST light curves for Luhman 16 
binned per 10 min intervals. For each light curve, the best-fit 
global model (see Sec. 3) is over-imposed in red. Gaps for night 
#1 and #5 correspond to cloudy conditions. The gap for night #2 
corresponds to the observation of another target. 



tational modulation was based on a division of the brown dwarf 
into 10 longitudinal slices. For each of them, the surface flux was 
assumed to be constant during each night, the measured flux be- 
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Fig. 2. Globally normalized TRAPPIST differential photometry for Luhman 16 (top) and for one of the comparison stars used in the 
reduction (bottom, shifted along the y-axis for the sake of clarity), unbinned (cyan) and binned per interval of 30 min (black). The 
standard deviation of the binned light curves are 2.2% and 0.1% for Luhman 16 and the comparison star, respectively. 
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Fig. 3. Differential photometry for the second part of night #10, 
binned per 10 min intervals, obtained with an aperture en- 
compassing both components of the binary (filled green cir- 
cles), and with apertures encompassing only the PSF center 
of Luhman 16A (black empty symbols) and Luhman 16B (red 
vertical bars). The amplitude of the variations is amplified on 
Luhman 16B, indicating that they originate from this T-dwarf. 



ing modeled by a semi-sinusoidal function falling to zero when 
the slice disappears from view. To this rotational model, we 
added a baseline model accounting for (i) the flip of the equato- 
rial mount of the telescope at the meridian, putting the stellar im- 
ages on different pixels and thus possibly creating small offsets 
in the differential photometry, (ii) a 2 -order airmass polyno- 
mial aiming to model the differential extinction curvature due to 
the much redder color of the target compared to the comparison 
stars, and (iii) a A th -order time polynomial representing the low- 
frequency variability of the system, including the evolution of 
the patterns from one rotation to the other. In this global model, 
the only two perturbed parameters in the MCMC were the rota- 
tion period and an arbitrary phase, the solution for the remaining 
parameters being obtained by linear regression at each step of 
the Markov chains (see G12 for details). Two MCMC chains of 
100,000 steps were performed to probe efficiently the posterior 
distribution of the rotational period. 



Our MCMC analysis gives P rot - 4.87±0.01h, producing an 
excellent fit between the model and the data (see Fig. 1). While 
a periodogram of the residuals reveals no additional period, we 
performed an additional analysis by adding a second rotational 
model to the MCMC. This analysis also failed to identify a sec- 
ond period. The Bayesian Information Criterion (Schwarz 1978) 
significantly increased by +360, indicating that a model includ- 
ing a single rotation period is a more likely representation of 
the data. We thus conclude that only one of the two BDs domi- 
nates the photometric variability since a same rotational period 
for both BDs is unlikely. 

A new photometric reduction of a fraction of our images us- 
ing the IRAF/DA0PH0T ALLSTAR PSF-fitting software (Stetson 
1987) revealed that the B component is -0.1 ±0.1 mag brighter 
than the A component 

We then attempted to determine the origin of the signal. We 
extracted the fluxes of both components by aperture photometry, 
fixing the aperture centers on the known positions of the two 
BDs. The aperture sizes were chosen to encompass only one PSF 
center. Using only data showing both a large signal and having 
the smallest and most stable PSF widths, we find that a larger 
amplitude is visible on the T-dwarf, while no significant signal 
is obtained for the L-dwarf. This is illustrated for a fraction of 
night #10 in Fig. 3. From these results, we conclude that the 
detected quasi-periodic variability originates from the T-dwarf 
Luhman 16B. 



4. Discussion 

From their new spectroscopy and from the relationship of 
Stephens et al. (2009), K13 attribute to Luhman 16A and B com- 
ponents effective temperatures of 1350±120K and 1220±110K, 
respectively. Such very low temperatures make atmospheric con- 
densates the most likely source of the observed rotational vari- 
ability (see discussion by Radigan et al. 2012 for the T1.5 BD 
2M2139). The most surprising feature of the patterns reported 
here is their fast evolution from night to night. To our knowledge, 
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this is the first time that such rapid evolution is reported for the 
cloud coverage of an L/T transition BD, making Luhman 16B 
an actual 'Rosetta Stone' for the study of BDs atmospheres. 
Future multi-band and high-cadence spectroscopic time-series 
will be able to observe while clouds condensate and are then dis- 
persed, giving unprecedented access into the physico-chemical 
processes acting in such atmospheres. 

Considering our measured Amag = 0.1 ±0.1 between both 
BDs, the variability amplitudes visible in Fig. 1 are thus diluted 
by a factor ~2 by the blend of both PSFs in the TRAPPIST 
images. The actual maximal peak-to-peak variability amplitude 
of Luhman 16B should thus be >20%. Similar amplitudes were 
observed around Ifim for at least another early T-dwarf, the 
T1.5-type BD 2M2139 (Radigan et al. 2012). For the L-dwarf 
Luhman 16 A, we do not detect any variability, but our sensitivity 
is limited by the partial resolution of both components in our im- 
ages. High-cadence photometry with a better spatial resolution 
will be need to thoroughly assess its photometric variability. 

Binary stars have been used from the dawn of modern astron- 
omy as a means to control some variables, while leaving others 
free. In our case both BDs probably formed together out of the 
same fragment; they should have their age and composition on a 
par. What differentiates them is their mass. This affects their re- 
spective effective temperatures and densities. Although only 100 
K apart, that one shows clouds breaking while the other does not 
is an indication of how sharp the boundary from a stable to an 
unstable atmosphere is, and of how closely tuned it must be to 
effective temperature and scale height. 

Field BDs are very interesting targets for exoplanets tran- 
sit searches (Blake et al. 2008; Bolmont et al. 2011; Belu et al. 
2013), as their small sizes make possible the detection of terres- 
trial planets with photometric precisions at the mmag level sim- 
ilar to the ones reported here for Luhman 16. A careful visual 
inspection of the residuals light curves and a search for periodic 
box-like patterns with the BLS algorithm (Kovacs et al. 2002) 
failed to detect any transit-like signal. Despite the blend of both 
BDs, our sensitivity is good enough to partially probe the terres- 
trial regime. Fig. 4 (top panel) shows the residuals of our global 
modeling compared to the expected transit depths for a 2 R @ ra- 
dius planet, for both BDs, assuming for each of them a Q.IR Q ra- 
dius. The corresponding transit is firmly discarded by our data. 
For orbital periods smaller than the mean duration of our runs 
(~9.5h), our detection threshold goes down to Earth-size planets, 
as can be see in Fig. 4 (bottom panel). This shows that intensive 
campaigns like the one described here targeting the most nearby 
field BDs with small to medium-sized ground-based telescopes 
or, better, with an infrared space facility like Spitzer (Triaud et 
al. in prep.) could efficiently assess the frequency of close-in ter- 
restrial planets around BDs, possibly detecting Earth-sized plan- 
ets amenable for atmospheric characterization with, e.g., JWST 
(Belu etal. 2013). 
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Fig. 4. Residuals of our global modeling, unfolded and binned 
per 10 min intervals (top), folded with P - 8h and binned per 
5 min intervals (bottom). For the unfolded residuals is shown as 
a red horizontal line the expected transit depth for a 2R @ planet 
orbiting one of the two BD, assuming a 0.1R Q radius and the 
same brightness for both BDs. The same is done for the folded 
residuals, but assuming a IR @ planet. 
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